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I. INTRODUCTION 

Entanglement, as a term of joint quantum coherence, 
is one of the most intriguing elements of quantum me- 
chanics and it is crucial in quantum information tasks 
[l[ . However the existence of an interacting reservoir or 
environment that leads to decoherence and/or disentan- 
glement places an obstacle to the maintenance of 
joint quantum coherence during any dynamical process 
|5| . Thus the study and control of entanglement dynam- 
ics has received wide attention in recent years (see re- 
views )■ There have been studies of entanglement 
dynamics from many points of view. Examples involve 
open system treatments [9rJl3j or closed quantum sce- 
narios such as cavity QED systems [T3 - |2j1 | , spin systems 

dum, etc - 

Many interesting and sometimes surprising findings 
such as entanglement sudden death d, [l(| El, sudden 
birth [12], revivals [l3l [Til l 20 | . dynamical relations with 
quantum state transfer 22-2J] , and other exotic types of 
entanglement evolution have been reported. Such inter- 
esting phenomena accompany the idea of tracking entan- 
glement as a carrier of quantum information [lll426l |. a 
generalization of entanglement swapping 30] . 

One consequence has been the discovery of examples 
of non-trivial "information conservation" among three 
or more parties [H, [III [26|, arising in cases of suffi- 
ciently symmetric interaction Hamiltonians, or special 
initial states, or reservoirs that are sufficiently small that 
their state evolution can be followed in detail, such as in 
perfect-mirror closed-system cavity QED [l7H2l| and in 
spin systems (23rl2fj| . However true reservoirs are com- 
plex and difficult to follow, especially if mixed state con- 
siderations are important. No general closed-form rules 
of entanglement transfer are known in such cases. 

In this paper we revisit quantum information flow from 
a different perspective and derive a new class of entan- 
glement constants of motion. Our approach employs am- 
plitude channel dynamics and avoids information loss by 
tracing, while remaining open to non-Markovian as well 



as Markovian reservoir behavior. We note that the sys- 
tem of experimental interest, which may be one or more 
qubits, is almost always prepared in a pure state if possi- 
ble, and frequently the method of preparation produces 
a pure entangled state. This means that the system it- 
self is in a mixed rather than a pure state. We assume 
that the entanglement during state preparation, causing 
the mixedness, arose via interactions that have ceased 
prior to the beginning of a period of interest at t = 0. 
This period of interest could simply be intended for quan- 
tum memory preservation or for specific state manipula- 
tions. The static disengaged nature of the prior entan- 
glement partner, and also its lack of specificity in our 
treatment, reduce it to a vague background object in any 
further qubit evolution, and for this reason we label it the 
"Moon" . 

A general sketch of our scenario is given in Fig. [T] 
Unit A is taken as a two-level system (qubit) and unit 
a as a separate quantum system of arbitrary dimension 
interacting with it, nominally a reservoir. The Moon 
M , i.e., the non-interacting, unspecified, and completely 
static background, is entangled via an earlier preparation 
stage with A. 

There obviously remains a wide choice for systems act- 




FIG. 1: A general sketch of our scenario. The bubble circles 
the system of interest A and leaves everything else out. The 
dashed line indicates its entanglement, but not interaction, 
with the unspecified background "Moon" M. The arrow rep- 
resents interaction between A and an arbitrarily-dimensioned 
unit a, which can be the quantum vacuum reservoir, a single 
mode cavity, an XY spin chain, etc. 
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ing as environments that promote evolution of the sys- 
tem of interest after t — 0. We will illustrate a range 
of possibilities with concrete results in various specific 
interaction contexts: spontaneous emission |3lj . Jaynes- 
Cummings ( JC) cavity dynamics [32[ , and XY spin chain 
interactions [33[. These present very different physical 
situations and interaction mechanisms, and lead to dis- 
tinct entanglement dynamics, but they all react similarly 
to the initial Moon entanglement. Our linked informa- 
tion constants arise from amplitude channel dynamics 
but do not rely on symmetries of the Hamiltonian or of 
any special initial state, in contrast to the cases in some 
previous work [26[ . 



II. SCHMIDT ANALYSIS OF ENTANGLEMENT 

In this section we address our approach to entangled 
state analysis. The Hamiltonian of our scheme reads 



H = H A + H a + H Aa + H 



(1) 



where Ha, H a and Hm are the Hamiltonians of the 
qubit system A, "reservoir" unit a and the previously- 
interacting Moon M respectively; and HAa denotes the 
only existing interaction, that between A and a. 

We start from the A-M entangled preparation state, 
i.e., the joint superposition state 



|^am(0)) = cos0|e)|mi) + sin%)|m 2 ), 



(2) 



where |e) and \g) are the excited and ground states of 
our qubit system A, and 



| mi) = ^2<ut\Mi), and 



m 2 ) 



(3) 



are two normalized Moon states with {|M/)} defined as 
a complete basis set for M. It need not be the case 
generally, but we assume in our example that the Moon 
states \m\) and \1TI2) are orthogonal, and because M is 
not interacting with either A or a they are effectively 
static: 



A-a part can be separated from M. Therefore we will 
only need to focus on the A-a dynamics when we study 
the time evolution. 

Before we proceed to the time dependent state in var- 
ious specific models in the following sections, we will 
first discuss the initial Moon entanglement. As we know, 
the pure-state relation between the two sides of any bi- 
partition is an R x S dimensional matrix (where R and S 
can be any numbers or infinite) that may connote entan- 
glement, but in any event permits a Schmidt-type decom- 
position of the joint state [33 - l36j . We use the Schmidt 
parameter K introduced by Grobe, et al., [34| as our 
quantitative measure of entanglement, where K is not 
simply the dimension of the space [l| but rather relates 
to the number of Schmidt modes that make a significant 
contribution to the state. Therefore we name this param- 
eter K the "Schmidt weight" from now on. The range 
of this Schmidt weight, N > K > 1 (JV is the effective 
dimension of the space) corresponds to the concurrence 
range 1 > C > 0, when concurrence [37j is also appli- 
cable. The upper and lower ends of both ranges denote 
maximal and zero entanglement, respectively. 

The Schmidt weight between two parties a and (3 of a 
general pure state \ip a p) is defined as 



K 



k 



(6) 



where these A^s are the non-zero eigenvalues of the re- 
duced density matrix for either system, p a or pp [36j : 



Pli = Tr Q [p] = Tr Q \ip a p)(lp a p\ 



= cc*. 



(7) 



Here C is the coefficient matrix connecting the two sep- 
arate arbitrary complete bases \n) and \p) of systems a 
and /3 respectively, with 



\lpap) = 2J C ( n > M)l«a) ® IM- 



(8) 



The square roots of A^s are also the coefficients of the 
usual Schmidt decomposition [lj 

(9) 



(m 2 (i)|mi(t)) = (m 2 |mi) = 0. 



(4) 



We adopt a conventional approach to the interacting 
partner system labelled a, assuming that it is separa- 
ble from the qubit A at t = 0, just as in the conventional 
treatment of a reservoir in quantum open system dynam- 
ics ■ Therefore the entire initial state can be written 
as 

\ipAaM(0)) = (cos0|e)|mi) +sin%)|TO 2 )) <g> |0i), (5) 

where \<fii) = J2k ik\k) is a normalized state of unit a, and 
{| A;)} is a complete basis for a. Usually \<j>i) is the ground 
state of part a. We note that since M is not interacting, 
the evolution of the states |e)|^i) and |g)|<£i) are driven 
only by the Hamiltonian Hau, i-e., the dynamics of the 



where \ff?) and \g%) are the orthonormal Schmidt states 
satisfying = (fff | 5 f,) = 6 kk ,. 

Since a general pure state is usually in some arbitrary 
basis other than the Schmidt basis, it is natural for us to 
follow the coefficient matrix procedure to calculate the 
Schmidt weight ([6]) . Accordingly we note from the initial 
state (JSJ) that the coefficient matrix for the Moon M in 
the basis of \mi), |m 2 ), (ma), and the interacting 
partner a in the basis {|fc)}, is an 00 x 00 matrix which 
is given as 



C 



M 



/ i\ cos( 


V : 



Ik cos f 






i\ sin( 





ik sin 9 



(10) 
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where the two-dot sign represents the elements 
ik cos9 for all the other fcs, while the three-dot sign 
represents empty rows or columns of zeros for the infi- 
nite number of remaining matrix elements. The reduced 
density matrix however is simply 



Pm — CmC m = 



cos 





sm 



(ii) 



We note that the qubit has reduced the effective inter- 
action space of the Moon to a two dimensional subspace, 
which means that in this context a two-state subspace of 
the Moon is in fact quite general. The non-zero eigen- 
values of the above matrix are obvious, and the resulting 
Schmidt weight Km denoting the entanglement between 
the Moon the remainder is given as 



K 



M 



cos 4 9 + sin 4 9 



(12) 



Since the Moon M is not interacting, its internal dy- 
namics only amount to a local unitary transformation 
38], which will not affect the entanglement between M 
and the rest, so Km is independent of time. We note that 
there is no Moon entanglement (Km = 1) when 9 = it/2 
or 0. In these cases the initial state is a trivial product 
state. Otherwise Km > 1- It is particularly interesting 
when the Moon restricts the A-a dynamics and acts as a 
monitor of the entire entanglement flow. The following 
sections take a few specific examples of the A-a inter- 
action to show the role of Moon entanglement in their 
particular entanglement information dynamics. 



III. SPONTANEOUS EMISSION 

In this case the qubit system A is a two level atom 
and unit a is the quantum vacuum reservoir consisting of 
the continuum of photon modes. The atom will of course 
decay to its ground state asymptotically and irreversibly, 
while one photon is emitted [3l[ HH . We write the Hamil- 
tonian in the usual way as a sum of atom and reservoir 
contributions: 



H A = ^huJA<J z A 



and H„ 



y^^huikajak- 



(13) 



Here a\ is the usual Pauli matrix, and the usual bo- 
son operators represent the reservoir with a continuum 
of modes, where a\ and at denote the standard cre- 
ation and annihilation operators respectively, and uja 
and LUk are the atom and reservoir frequencies. Here 
k = 1, 2, 3, oo, labels the infinitely many modes. The 
interaction Hamiltonian is also standard: 



H, 



E 

k 



H9*k<7 A a k+ 9k^a k ), 



(14) 



where crjj, er^ are the usual raising and lowering Pauli 
operators for the two level system A, and the gkS are 
coupling constants between the reservoir and the atom, 
for which fundamental expressions are well known (39j | 



\9k\ 



2he V 



rd? n COS 2 9. 



(15) 



Here 9 is the angle between the atomic dipole moment 
dio and the electric field polarization vector e^, and V 
is the quantization volume. According to our generic 
description in Eq. ([5]), the initial state can be rewritten 
in the spontaneous emission case as 



IV'AoAf(O)) = ( cos0|e)|mi 



3 in%)|m 2 )J ®|0>, (16) 

where |e) and \g) are the excited and ground state of the 
two level atom, and we have defined \4>i) = |0), indicating 
that all the reservoir modes are in their vacuum states. 
With the help of the Weisskopf-Wigner treatment [3l"l.l39| 
we will find 

\ipAaM{t)) = cos9\mi(t))(c e (t)\e)\0) + J2 c k(t)\g)\lk)) 

k 

+ sm9\m 2 (t))\g)\0), (17) 



where the coefficient c e (i) 



' r At/2 with 



d 2 w^/37reo?ic 3 as the natural line width, denotes 
that there is one photon in the reservoir mode fc while all 
the rest of the modes are empty and the coefficients 



ck(t) = gk 



1 



(18) 



are time dependent. The probability to find the atom 
in the excited state is P e {t) — \c e (t)\ 2 = e~ TAt , which 
decays to zero asymptotically and irreversibly. 

With the dynamical state (fTTj) we can begin to cal- 
culate the Schmidt weight KA(t), or K a (t), representing 
the entanglement between qubit A, or vacuum reservoir 
a, and their corresponding remainders. As defined in the 
last section, the coefficient matrix between A and the 
remainder for the time dependent state (|17p is given as 



C A 



c e (t) cos( 







... 

sin 9 a(t)cos9 .. Ck(t) cos 9 . 

(19) 

where the two-dot sign represents Ck(t)cos9 for all 
the other fcs and the three-dot sign again represents 
empty rows or columns of zeros. Then the reduced den- 
sity matrix is simply a 2 x 2 form 

' \c e (t) | 2 COS 2 9 

J2k Mt)| 2 cos 2 <9 + sin 2 ( 

. (20) 

Now according to the definition in Eq. (|6]), we imme- 
diately have the qubit entanglement Xa(*) given by the 
expression 



PA 



K A (t) 



[2 cos 2 9e 



-r A t 



(21) 
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We note that at t = 0, 



K A (0) 



cos 4 9 + sin 4 1 



(22) 



a finite number that is naturally the same as the con- 
stant Moon entanglement (|12p . As time goes on, the 
probability of the atom in the excited state P e (t) decays 
gradually, and at t = oo the probability is completely 
transferred to the ground state and leaves the atom in 
a product state with its remainder system, which means 
eventually A is disentangled from the rest of the universe, 
and the Schmidt weight Ka(oo) = 1. Fig. [2] illustrates 
the behavior .Fd(£) as a function of t at four different 9 
values. We note that in the region when sin 2 9 < cos 2 9, 
Ka(£) starts from a finite value, evolves to a local max- 
imum and then decays irreversibly to 1 as is shown in 
Fig. [21(a) and (b). However when sin 2 9 > cos 2 9 as is 
shown in Fig. [2] (c) and (d), KaII) decays directly and 
irreversibly to 1 . 

Now let us focus on the reservoir entanglement. From 
the time dependent state (fTT)) we see that the coefficient 
matrix of reservoir a is given as 



/ sin6» c e (t)cos9 ... \ 
Cl {t)cos9 



C a = 



Ck(t) cos 9 



V 



(23) 



/ 



Then the reduced density matrix is given by p a = C a Cl 
with an infinite number of non-zero eigenvalues 



A 



|c e (i)| 2 cos 2 ( 
|c fc (t)| 2 cos 2 i 



sm 



(24) 
(25) 



for k = 1, 2, 3, oo. Now from the definition ([6]) we find 
the reservoir Schmidt weight: 



K a {t) 



[2cos 2 6»(l-e- r ^)-l] 2 + l' 



(26) 



Obviously the reservoir is initially not entangled. Then 
its entanglement gradually increases, and at time t = 
(\n2)/T A we find K A (t) = K a (t). When time goes to 
infinity we note 



K a {co) = K A (0) = 



1 



cos 4 9 + sin 



(27) 



That is, the final reservoir entanglement K a {oo) equals 
the initial qubit entanglement Ka(0). Fig. [2] plots the 
behavior of K a (t) as a function of t for various 9 val- 
ues. We note that in the region when sin 2 9 < cos 2 9 as 
is shown in Fig. [5] (a) and (b), K a {t) starts from zero 
entanglement, reaches a maximum and then evolves to a 
finite value Ka(0) in the end. In the opposite region of 9 
as shown in Fig. [5] (c) and (d), it increases directly and 
irreversibly to the value Ka(0). 




FIG. 2: Time dependence of Schmidt weights KA(t) and 
K a (t) for the spontaneous emission dynamics at different val- 
ues of 9. In each of the four plots, the blue solid and the red 
dotted lines denote .fd(t) and K a (t) respectively, and the x 
axis represents time t with unit 1/Ta. Complementary be- 
havior of the two Schmidt weights is shown clearly in plots 
(c) and (d) where sin 2 6 > cos 2 6. 



To summarize, the qubit entanglement K A (t) starts 
at a finite value and decays completely to zero entan- 
glement, while K a (t) starts from no entanglement and 
eventually inherits the exact amount of the qubit 's ini- 
tial entanglement K A (0). This is exactly equal to the 
Moon entanglement Km, so one can see that the un- 
known Moon's entanglement (|12p has jumped into the 
picture. It is constant itself, but it acts as a kind of 
buffer to restrict information flow to and from A. In an- 
other way of speaking, we could say that there is only 
a certain amount of "free" entanglement able to be ex- 
changed, which is determined by the Moon. 

To take a further step and without loss of generality, 
we now assume sin 2 9 > cos 2 9 for convenience. Then we 
find two equalities connecting each of KA(t) and K a (t) 
to K M : 



K A {t) 



-T A t 



At )cos'9, (28) 



K a {t) 



2 
Km~ 



I = 2e 



- TAt cos 2 1 



(29) 



We note that both K A (t) and K a (t) are controlled by 
the Moon in a non-linear way. This control leads to a 
novel conservation relation between A and a in a pairwisc 
fashion: 



K A {t) 




K a (t) 



1 = 1 



~7? 1- (30) 
-Km 



Although K A if) and K a (t) are time dependent quantities 
they combine in this way to a constant determined only 
by the Moon entanglement Km- 
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Because of the restriction of the Moon entanglement 
we note from the conservation relation ([50)1 that in this 
region sin 2 9 > cos 2 9 or 9 g [tt/4, 37r/4], the decreasing 
of K A (t) is accompanied by the increasing of K a (t) as 
is shown in Fig. [2] (c) and (d). To see quantitatively 
this complementary relation let us take 6 = 7r/4 as an 
example. Then the time dependent qubit and reservoir 
Schmidt weights simplify to 



K A (t) 
K a (t) 



(1-e 
2 



-r A t\2 



,-2T A t i i • 



(31) 
(32) 



Obviously these two equations for K^it) and K a (t) de- 
pend on time in an opposite way. It is interesting to note 
that the parameter Ta, which represents a collective cou- 
pling between the atom and the reservoir modes, is also 
controlling the two entanglements inversely. This is be- 
cause of the Moon entanglement Km, which acts as a 
buffer to both entanglements KA(t) and K a (t) but sub- 
stantially in a opposite way through Ta (see Eqs. (|2"5)l 
and (j29|) ). 

We remark that for sin 2 9 < cos 2 9, a similar relation 
to Eq. ([30| can be achieved with only a modification 
of signs. In this case KA(t) and K a (t) are not always 
complementary any more as Fig. [2] (a) and (b) show 
for wide regions of t. However KA(t) and K a (t) still 
stand in a time-invariant relation similar to fl3Q} and are 
connected only by Km- 



IV. JAYNES-CUMMINGS INTERACTION 

Spontaneous emission is an example of an irreversible 
process. In this section we will turn to a simple example 
when the A-a interaction is reversible, following the JC 
model [32]. Thus qubit A is still a two level atom while 
unit a is now simply a single mode lossless cavity. Local 
entanglement dynamics between the atom and the field 
in the JC model was first studied by Phoenix and Knight 
[lij by expressing the entangled atom-field state in terms 
of the eigenstates and eigenvalues of both the field and 
atomic operators, and revival physics (ioMi^ played a 
key role in the dynamics. Later, it was shown by Son, et 
al., fioj j that entanglement between two non- interacting 
qubits can be generated through the qubits' local inter- 
actions with their corresponding JC cavities that are ini- 
tially in an entangled two-mode squeezed state. Lee, 
et al., [13] showed that there is actually entanglement 
reciprocation between the two qubits and their corre- 
sponding continuous-variable systems such as JC cavi- 
ties. Recently, the JC model was revisited by Yonac, et 
al., [13, HH, ISO]; an d by Sainz and Bjork [l9j . to illus- 
trate the entanglement sudden death phenomenon Q , as 
well as to track the entanglement flow, and conservation 
relations were found by both groups [HI, HH . 

Here we will continue to track the entanglement in- 
formation in the JC dynamics, but in addition will ac- 



count quantitatively for the role of the non-interacting 
unknown Moon. The JC Hamiltonian is given as 

HAa = -fwJA^A + %(aV^ + aa\) + Tibj a a)a, (33) 

where a A , a\ are the usual Pauli matrices describing the 
two level atom A, while and a denote the standard 
creation and annihilation operators for the single mode 
cavity. The atom and cavity frequencies are loa and ui a , 
respectively, and g is the coupling constant between the 
atom and the cavity. For convenience we take the reso- 
nant condition when uja = uj a . 

Now from the generic expression ([5]) the initial state 
for the JC model can be written as 

\ipAaM (0)) = (cos0|e)|mi)+sin%)|m 2 )) ® |0>, (34) 

where |e) and \g) are the excited and ground state of the 
two level atom, and we have defined \<j>i) = |0) as the zero 
photon state of the cavity. From the Jaynes-Cummings 
treatment [32j we will have the time dependent state as 

\i>AaM{t)) = co S 9\m 1 (t))(e l "* t / 2h cosgt\e)\0) 

- ze-^ 2h smgt\g)\l)) 

+ S in0|m 2 (t))| 5 >|O), (35) 



where |1) means that there is one photon in the cavity. 
If we follow the same Schmidt calculations as in the last 
section, we will have entanglement KA(t) between atom 
A and the rest of the universe as 



K A {t) 



[2 cos 2 9 cos 2 gt - 1] 2 + 1' 



(36) 



That is, we find expression ([2~Tj) again, except that e~ TAt 
has been replaced by cos 2 gt. This is just the replacement 
of one formula for excited state probability by another, 
as the nature of the amplitude decay channel requires. 

While the initial qubit entanglement Ka(0) again has 
a value equal to the Moon entanglement ([12"]) , in the JC 
dynamics K a (t) has a period of r = 7r / g instead of decay- 
ing irreversibly as in the spontaneous emission case. We 
see at the half period time the atom loses all of its entan- 
glement: KA{t = t/2) = 1. Then it evolves to the initial 
value Ka(0) at t = r. Fig. [3]shows this periodic behavior 
of -fCt(i) plotted as a function of t at different 9 values. 
Recovery of atom-field and atom-atom entanglement in 
the JC dynamics was already shown previously in Refs. 
[14H20I ]. However, here our result shows a different type 
of entanglement recovery, because KA(t) denotes another 
type of entanglement, this time including the unspecified 
non-interacting Moon as well as the cavity. 

The cavity entanglement K a (t), 



K a (t) 



[2 cos 2 9 sin gt - l] 2 + 1 



(37) 



is also predictable if we look to ([26| and see that 1 



-r A t 



should be converted to sin gt because both are 
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FIG. 3: Time dependence of Schmidt weights KA(t) and 
K a (t) for the JC dynamics at different values of 9. In each 
of the four plots, the blue solid and the red dotted lines de- 
note KA{i) and K a (t) respectively, and the x axis represents 
time t with unit 1/g. Complementary behavior of the two 
Schmidt weights is shown clearly in plots (c) and (d) where 
sin 2 > cos 2 e. 



expressions for the ground state probability. We note 
that K a (t) is also periodic. It is initially not entangled 
with its remainder (K a (Q) = 1), and then increases with 
time. At t = ir/Ag, we have K a (t) = KA{t), and at the 
half period time we see that 



K a (t = r/2) - K A (0) = K M , 



(38) 



exactly the same as Ka(0) entanglement. Again Fig. [3] 
illustrates the periodic behavior of K a (t) as a function of 
t at various 9 values. When compared with the behav- 
ior of the qubit entanglement we see that the amount of 
entanglement Km has been completely transferred from 
KA(t) to K a (t) at the half period time t = t/2. After 
this, however, the entanglement is repeatedly transferred 
back and forth between K A (t) and K a (t). This is the ma- 
jor difference from the spontaneous emission case where 
the reversible process is absent. 

Again we work in the sector when sin 2 9 > cos 2 9 
for convenience and see that both of the entangle- 
ments K A (t) and K a (t) are restricted by the constant 
Moon entanglement Km in the following non-linear time- 
dependent way: 



^-1-^-1 = 2sinVcos 2 *, (39) 



_^-l-^J--l = 2cosVcos 2 .. (40) 



This periodic time dependent control of the two Schmidt 
weights by the Moon entanglement is different from the 
spontaneous emission case. However, the two equalities 
also lead to the same generic entanglement conservation 



A 



N-l 



-O 

N 



J, 



FIG. 4: Scheme of qubit A interacting with an XY chain. The 
big blue dot represents the qubit A while the N small dots 
represent the spins in the XY chain. The arrow indicates 
interaction between A and the first spin in the chain with 
coupling constant J a while the solid segments connect the 
nearest neighbor sites in the chain with coupling constant J. 



relation 



K A {t) 




K a {t) 



1 = 1 



2 
K M 



(41) 



Therefore in the JC model case the time dependent 
Schmidt weights KA(t) and K a (t) are also restricted by 
the constant Moon entanglement Km- We see clearly 
here that the decrease of K A (t) is accompanied by the 
increase of K a (t) and vice versa as is shown in Fig. [3] (c) 
and (d). To show quantitatively we again take 9 = it /A 
as in Fig. [3] (c) to follow this complementary relation of 
the two entanglements: 



K A (t) 
K a (t) 



sin gt + 1 
2 



(42) 
(43) 



cos 4 gt + 1 ' 

which are the exact analogs of (f3T|) and (j3"2")l . 

V. XY SPIN INTERACTION 

We now move to a condensed matter context and take 
a final example when the A-a connection is a Heisenberg 
exchange interaction or spin-spin interaction. Here qubit 
system A is a spin one-half particle while unit a is now 
an TV-spin XY chain [33| (see Fig. 2]), a simplified model 
for strongly correlated materials such as ferromagnets, 
antiferromagnets, etc. 

The first studies of entanglement flow in spin chains fo- 
cused on few-qubit chains (N < 6) and the W state [27l ]. 
and also entanglement dispersion in long chains (N 3> 1) 
[28j |. Amico, et al. [29[ studied the propagation of a 
pairwisc entangled state through an XY spin chain, and 
found that singlet-like states are transmitted with higher 
fidelity than other maximally entangled states. Here we 
also focus on the entanglement dynamics, not to trans- 
port the entanglement, but to track the information flow 
by taking into account the role of the entangled Moon. 
The interaction Hamiltonian of our scheme is given as 



H Aa = Ja{?\vx +v a (J i) 



N-l 



^2 J ( (7 i <T n+l + a n O'n+l)) 



(44) 



where a ± are the usual Pauli matrices describing the 
spins, J a is the coupling constant between spin A and 
the first spin a± of the XY chain, and J is the coupling 
constant between the nearest neighbor sites inside the XY 
chain. Now we take J a — J for convenience. This model 
can be transformed through a Jordan- Wigner transfor- 
mation |43l| into a set of free fermions (see for example 
Ref. [44]) and thus can be solved exactly (33[. From 
the perspective of the Jordan- Wigner transformation, the 
XY model is equivalent here to a free fermion hopping 
model or Tight-Binding model describing phonon sys- 
tems. 

For the XY Hamiltonian Hao, the exact N + 1 eigen- 
states are given as 



\k) 



N ■ 



A"7T 



It) |0> 



N 



N ■ 



N + 2 

(n + 1) kir 



N + 2 



II) I In), (45) 



where |f) and ||) are the spin up and down states for 
our qubit A, and |l n ), |0) are the states of the XY spin 
chain with |1„) indicating there is a spin up at the nth 
site while all the rest are in the spin down state and 
|0) meaning all the sites from site 1 to site N are in the 
down state. Here k = 1,2, 3, N+l represent the iV + 1 
eigenstates. The corresponding eigenvalues are 



Ek = 2 J cos 



k'TT 



Then the evolution operator can be written as 

JV+l 



U Aa {t) = e 



-iE k t 



\k)(k\. 



(46) 



(47) 



fc=i 



Again from the generic initial state ((5|) we have here 
for the XY model 

\^Aa M (0)) = ( cos0 |t) |mi) + sin0|J.) \m 2 )) <8|0), (48) 



where we have defined \4>i) = |0) to represent all the N 
spins in the XY chain that are in the down state. Then 
the time dependent state can be achieved as 



|^AoAf(*)) 



cos 



71=1 



+ sin0|ma(t))| |)|0). 
where we have defined 

N+l 



0|m 1 (O)( Ce Wlt)|O) + ^c„(i)ll)|ln 

(49) 



, s 2e' lEkt 



k=l 



N - 



(V4-1 



kir 
N + 2 



N + 2 



(50) 



Since c e (t) and c n (t) are complicated expressions for 
arbitrary number N, here we take N = 10 as an example 
to illustrate their properties. Then we have 

c e (i) = — 2 + 3cos(Jt) + 2cos(V2J£) + cos(v / 3Jt) 



+ (2 + ^)008 

+ (2 - VS) cos 



(V3-l)Jt 

72 
(V3 + l)Jt 



V2 



(52) 



We note that the five cosine functions have five differ- 
ent periods and the ratio of any two periods is irrational. 
Therefore the five quantities will not have a common pe- 
riod, which means that c e (t) will oscillate all the time but 
without a fixed period. Now we define f(J,t) = \c e (t)\ 2 
and note that it can vary from to 1 . There are infinitely 
many solutions for /(J, t) = as a function of time t, say 
t = Ti, with i = 1, 2, 3, oo. 

If we follow the same Schmidt calculations as in the 
last two sections we will find the Schmidt weight KA(t) 
between the qubit spin A and the remainder as 



K A (t) = 



[2f(J,t) cos 2 9- lp 



1 



(53) 



We note that the qubit entanglement KA(t) is also os- 
cillating as determined by f(J,t). As the amplitude 
channel requires, it starts at the familiar same value 
Ka(0) = Km, and in this example evolves to zero en- 
tanglement at the time points r^. After each of these 
zeros, KA(t) will increase to a local maximum point and 
then decay to 1 again at the next time point r^+i. Fig. [5] 
illustrates this particular behavior of Ka(J-) as a function 
of t at different values of 9. Such aperiodic behavior is 
intermediate to the previous two examples showing irre- 
versible decay and periodic oscillation, and is expected 
on the basis of the irrationally related spin-chain eigen- 
frcquencies. 

Now we come to the XY chain entanglement K a {t) 
representing the entanglement between the chain and its 
remainder, i.e., the end spin A and the Moon M. It is 
related to Ka in the usual way. We just replace f(J,t) 
by 1 — /(J, t) and obtain: 



K a {t) 



(54) 



[2cos 2 6»(l-/(J,t))-l] 2 + l' 

So the chain entanglement also oscillates with /(J, t). In 
general, the entanglement will be transferred back and 
forth between iGt(i) and K a (t) with Km the upper limit 
of entanglement that can be transferred just as the pre- 
vious two cases. 

As expected, restrictions on entanglement flow follow 
the previous examples. In the regime sin 2 9 > cos 2 9 we 
can simply repeat relations (|3"9l and PD|) by replacing 
cos 2 at with f(J,t): 
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FIG. 5: Time dependence of Schmidt weights KA(t) and 
K a (t) for the XY spin chain dynamics at different values of 
0. In each of the four plots, the solid and the dotted lines de- 
note KA(t) and K a (t) respectively, and the x axis represents 
time t with unit 1/J. Complementary behavior of the two 
Schmidt weights is shown clearly in plots (c) and (d) where 
sin 2 6 > cos 2 6. 

y / 5^1-^ 7 l = W, *)«»»*,. (56) 

Naturally, the same non-linear conservation relation (j4Tj) 
is recovered, and again the A and a entanglements behave 
complementarily, this time as a function of Jt as shown 
in Fig. M(c) and (d). 

VI. SUMMARY 

In summary we have studied entanglement information 
flow from the perspective of a dynamical qubit A in an 
initially mixed state, a state that was generated by an en- 
tanglement associated with a prior process, which we can 
loosely assign to an experimental preparation stage. Us- 
ing Schmidt-decomposition rather than master-equation 
analysis, we derived conservation statements for the sep- 
arate degrees of quantum entanglement of the qubit and 
of its interacting reservoir, and showed their relation to 
the entanglement of the unspecified background party 
we called the Moon, which was initially entangled but at 
t = ceased to interact with either the qubit A or its 
environment a. 

The new forms of entanglement conservation relations 
are nonlinear connections between quantum memories, 
dependent on the restrictions implied by amplitude flow 



channel dynamics. One can say that the channel's en- 
forcement of excitation number conservation in the qubit- 
reservoir interaction is the root cause of the entanglement 
and its flow. This is closely analogous to the continu- 
ous entanglement between transverse momenta in spon- 
taneous parametric down conversion, which arises from 
the enforcement of simultaneous momentum and energy 
conservation on the two-photon amplitude in the creation 
of the signal and idler photons. 

Although unspecified, and ignored in previous open 
system analyses, the Moon can be assigned responsibil- 
ity for the initial impurity of the qubit state. The three- 
part total universe (A + a + M) was bi-partitioned three 
ways in order to evaluate the respective Schmidt weights, 
as indicators of entanglements in three specific interac- 
tion models (spontaneous emission, JC interaction, and 
XY spin chain). These were analyzed to illustrate the 
flow of quantum information in different contexts, includ- 
ing both discrete and continuous versions of the reservoir 
system labelled a. Although the influences on individ- 
ual entanglements differ in various ways, the amplitude 
flow common to them produces entanglement conserva- 
tion relations in the same form. One can say that the 
non-specified Moon retains a kind of influence on the 
system of interest whether we are "looking" (through 
interaction) at it or not. The qubit can feel, through 
the entanglement conservation relation but not through 
interaction, that the Moon is there. 

There can be interesting consequences when the Moon 
also has a significant dynamical evolution, although still 
not interacting with A, because its entanglement with 
A can then be assigned to part rather than all of it. 
This discussion will be undertaken elsewhere (45J. Fi- 
nally we would like to comment on the inverse depen- 
dence of KA(t) and K a (t) on the interaction parameters 
as discussed at the end of our three examples. It will be 
particularly interesting if, for some systems, the interac- 
tion constant can be adjustable (e.g., the coupling con- 
stant of a spin-spin interaction) . Especially in the ther- 
modynamic limit interesting phenomena such as quan- 
tum phase transitions flfjl l47j may arise from changes of 
the interaction parameter. The behavior of the entan- 
glements in the vicinity of the critical point will be ex- 
tremely interesting (see for example [48] and references 
therein). 

We acknowledge helpful conversations with Profs. L. 
Davidovich and Ting Yu, and partial financial support 
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W911NF-09-1-0385, NSF PHY-0855701. 
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